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Q . Abstract 

"^ \ We develop a new method of implementing the Hartree-Fock calculations. A 

class of Gaussian bases is assumed, which includes the Kamimura-Gauss basis-set as 
^ . well as the set equivalent to the harmonic-oscillator basis-set. By using the Fourier 

rH \ transformation to calculate the interaction matrix elements, we can treat various 

interactions in a unified manner, including finite-range ones. The present method 
is numerically applied to the spherically-symmetric Hartree-Fock calculations for 
the oxygen isotopes with the Skyrme and the Gogny interactions, by adopting the 
C^ ' harmonic-oscillator, the Kamimura-Gauss and a hybrid basis-sets. The characters 

of the basis-sets are discussed. Adaptable to slowly decreasing density distribution, 
the Kamimura-Gauss set is suitable to describe unstable nuclei. A hybrid basis-set 
of the harmonic-oscillator and the Kamimura-Gauss ones is useful to accelerate the 
convergence, both for stable and unstable nuclei. 
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1 Introduction 

Since the invention of the secondary beam technology, numerous experimental data on the 
unstable nuclei have disclosed new aspects of the atomic nuclei. Remarkable examples are 
the presence of the nuclear halos and skins, and the dependence of magic numbers on the 
neutron (or proton) excess ^. It should be noticed that both are closely related to the 
properties of the single-particle (s.p.) orbits in the unstable nuclei. In order to understand 
these new phenomena, which have raised questions on some of our conventional picture 
of the nuclear structure, it is worthwhile reinvestigating the properties of the s.p. orbits 
in nuclei. 

Because the atomic nuclei are bound without an external field, a mean-field is necessary 
to obtain the s.p. orbits. The Hartree-Fock (HF) theory and its extensions, which are 
self-consistent approaches, will be a desirable tool to investigate the s.p. orbits from 
microscopic standpoints. In studying the structure of unstable nuclei, it is a practical 
problem how to treat numerically the wave- functions at relatively large r, because there 
may be a halo. Most of the methods employed so far are adapted to the nuclei with sharply 
decreasing densities at the surface. They are not necessarily eligible to reproduce the halo 
structure. Furthermore, an important physics problem lies in the effective interaction. Not 
many effective interactions have been used in the HF calculations of nuclei. The Skyrme 
interaction [Q has been popular in the HF studies, since the zero-range form is easy to 
be handled. A number of parameter-sets have been proposed for the Skyrme interaction. 
In the Skyrme interaction the non-locality in the nuclear interaction is approximated by 
the momentum dependence of the zero-range force. This approximation was justified by 
Negele and Vautherin via the density- matrix expansion |^. However, despite the success 
for the stable nuclei and its recent development 0], it has not been inspected sufficiently 
whether the first few terms of the density-matrix expansion give good description of 
the nuclei far from the /5-stability. In this regard, it is desired to deal also with finite- 
range interactions. The Gogny interaction [^ is the only finite-range interaction widely 
applied to the mean-field calculations. In almost all recent studies based on the Gogny 
interaction, the DIS parameter-set is employed. However, the DIS set has a problem 
which is revealed in the unstable nuclei [^. It could be important to consider various 
possibilities of the effective interactions in the mean-field calculations. 

In this article, we develop a new method to implement the HF calculations. The 
following two points will be kept in mind: (i) for the drip-line nuclei the wave-functions 
in the asymptotic region could be significant and therefore should be treated properly, 
and (ii) the method should have capability of handling various effective interactions, 
particularly some finite-range ones. Satisfying these two conditions, the method developed 
in this paper will be useful to study structure of unstable nuclei within the HF framework. 

2 Single-particle bases 

In the following discussions we assume that the s.p. orbits maintain the spherical sym- 
metry, for the sake of simplicity. The extension of the method to the symmetry-breaking 



cases will be straightforward. 

The HF calculations are implemented by solving the s.p. Schrodinger equation {i.e. 
the HF equation) iteratively. There are two well-known ways to solve the s.p. Schrodinger 
equation. One is to discretize the radial or spatial coordinate with a finite mesh, and to 
integrate the differential equation numerically. The other is to introduce a basis-set and 
to reduce the equation to an eigenvalue problem, by applying the matrix representation 
to the s.p. Hamiltonian. Unless the effective interaction has zero range, the HF equation 
becomes an integro-differential equation because the Fock term is non-local. This makes 
the mesh method to be cumbersome. On the contrary, we can store the two-body matrix 
elements in the basis method, and then the non-locality in the Fock term addresses no 
essential difficulty. Since we would deal with finite-range interactions as well as zero-range 
interactions, it will be advantageous to introduce a certain set of s.p. bases. 

Because dimensionality in practical calculations is necessarily finite, the calculated 
wave-functions more or less inherit characters of the original bases. Therefore the choice 
of the basis-set is important, and could depend on the system under discussion, in general. 
The harmonic-oscillator (HO) basis-set has been popular in describing the s.p. orbits of 
nuclei. However, while the HO set is indeed efficient in the stable nuclei, this is not the 
case for the drip-line nuclei, as will be shown in Section ^ The density of the drip-line 
nuclei slowly decreases for large r (radial coordinate). On account of the short-range 
character of the nuclear force, the asymptotic form of the density distribution should 
be exponential, e"^''*' [Q, where r] = ylME /Tic with the separation energy E. In the 
drip-line nuclei, the density in the asymptotic region could sizably contribute to physical 
quantities such as the rms radius and the binding energy. This is a sharp contrast to 
the stable nuclei. However, the exponential asymptotics is hardly expressed by the HO 
bases. It is noticed that the exponent depends on E, which is not obtained until the 
convergence in the HF calculation. We have to reproduce not only the exponential form 
but also the E dependence of the exponent, for the proper description of the asymptotics. 
It is commented that, in the mesh method, we need a large number of mesh points to 
reproduce the density distribution in the asymptotic region, as far as we keep the mesh 
size uniform. 

We consider the s.p. bases having the following form, 

¥.,,,„(r) = /2„,,(r)[rW(r)x.]^); 
Raij{r) = Ar„,,/+'^"exp[-(r/z/,)2]. (1) 

Here Y'^^^y) expresses the spherical harmonics and x<t the spin wave-function. We drop 
the isospin index without confusion. The index a indicates the extra power of r (pa), 
which is a non-negative integer, and the range of the Gaussian (//„), simultaneously. The 
constant Matj is determined as 



Af^e, = , , — — , (2) 

7ri./(2£ + 4p„ + l)!! V'^a 



so as for {(paejml'^aejm) to be unity. Since the bases of Eq. (|I]) are non-orthogonal between 
different a's, the Schrodinger equation leads to a generalized eigenvalue problem when 



these bases are applied. If we take pa = 0, 1, 2, ■ ■ ■ with a constant range z/^ = z/^ 



2h/Mu, the space spanned by these bases is equivalent to that comprised of the HO 
bases. Indeed, these bases coincide with the HO ones, whose radial part is given by the 
associated Laguerre polynomials of 2(r/z/^)^, if the Gram-Schmidt orthogonalization is 
carried out from the smaller p„ to the larger. Because all the HO bases have the common 
Gaussian factor e~^'^^'^"^ , the superposition of a limited number of the HO bases has 
the asymptotic form of e~^'^^'^'^^ again. This is the reason why the HO basis-set fails to 
describe the density (namely, the wave-function) asymptotics of the drip-line nuclei. On 
the other hand, Kamimura proposed a basis-set |^ in which z/^ is given by a geometric 
progression, while keeping pa = 0. This set of bases, which will be called Kamimura- Gauss 
(KG) basis-set in this article, has been shown to work efficiently in few-body systems ||10| , 



including loosely bound ones. Although each KG basis has the Gaussian asymptotics, 
the exponential decrease of the density at large r is appropriately described to a good 
approximation by the superposition of the Gaussians with various ranges. The present 
form of the bases (|I|) covers both the HO and KG bases. 

The form of Eq. ([l|) allows wider variety of basis-sets than the HO and the KG sets. 
An immediate possibility is a hybridization of the HO and KG bases. Another possibility 
may be stochastic selection of p^ and z/^, though it is not explored in this paper. 

The transformed harmonic-oscillator (THO) bases were developed to reproduce the 
exponential asymptotics in the density of the loosely bound nuclei, and were applied to 
the mean-field calculations with the Skyrme interaction |ll|] . In order to obtain the energy 



dependence of the exponent, the bases themselves are changed iteratively. It was shown 
that the THO basis-set gives an improved description of the properties of nuclei near the 



neutron drip-line, over the HO set. However, in Ref. |Tl| the actual calculation seems to 
depend on the characteristics of the zero-range interaction. It might not be easy to deal 
with finite-range interactions by the THO bases. 



The bases of Eq. (|l|) give the norm matrix of, for each {i,j] 

rim 



N^p — {<^a(.j'niWpijrn) 



{2i + 2p„ + 2pp + 1)!! (vpY'^-PP ( 2v^vp 



{2l + Ap^ + l)\\{2l + App + l)\\ ^^cJ \^l + ^l. 



I • (3) 



Under the £ and j conservation, the s.p. Hamiltonian matrix is given by 

^i? = {'^aljm\h\^mm) , (4) 

where h stands for the s.p. Hamiltonian. Suppose that \'ipntjm) is a solution of the s.p. 
Schrodinger equation, 

h\iJntjm) = ^ntj\'^nijm) ■ (5) 

By expanding \'ilJnijm) by the bases \'^aijm), 

IV'n^im) = Yl (^ni\V'^ijm) , (6) 



the Schrodinger equation (|) is represented as the generahzed eigenvalue problem, 

Since the norm matrix A^^^-^) is real symmetric, the Cholesky decomposition can be ap- 
plied, which is equivalent to the orthonormalization of the bases. Then the generalized 
eigenvalue problem is converted to the normal eigenvalue problem. 

When we deal with non-orthogonal bases, we have to be careful for the norm after the 
orthogonalization not to be too small. In particular, the bases of Eq. (|l|) could compose 
an over-complete set, as is obvious from the fact that the HO basis-set, which is equivalent 
to the set of the bases having a fixed z/q,, can already be complete. If the norm after the 
orthogonalization were too small, i.e. one of the bases were almost linearly-dependent 
on the other bases, a numerical instability could take place. This condition may pose a 
practical limit on the present method in choosing p^ and u^. 

3 Effective interaction 

The effective Hamiltonian for the nuclear mean-field theory consists of the kinetic energy 
and the effective interaction, 

H = K + V- K = Y.^, ^ = E%- (8) 

i i<j 

Here i and j are the indices of each nucleon. The s.p. matrix element of the kinetic term 
is calculated as 

{^aejm\-^\'f^ijm) 



1 



2Mp^up 



{21 + 2p„ + 2pp + 3)4^ - 2 {(£ + 2p,)^ + {£ + 2p^)^ 



+4 



2t + 2pa + 2p/3 + 1 2vaV(i 



^t$ ■ (9) 



It will be natural to assume the effective interaction Vij to be translationally invariant, 
except for the density dependence mentioned below. 

As stated in Section |l], we would consider various types of the two-body interaction. 
For the zero-range interaction like the Skyrme force, the s.p. Hamiltonian is represented 



in terms of the local densities and currents [y, 12|. It is fast to compute the matrix 
elements of the s.p. Hamiltonian via the local densities and currents. However, this is not 
the case for finite-range interactions. Hence we shall calculate the two-body interaction 
matrix elements and store them, as will be discussed in the subsequent section. 

The saturation must be fulfilled in the nuclear HF approach. This requires compo- 
nents other than the momentum- independent two-body terms in the central force |T3| . 



A density-dependent (or a three-body) interaction is usually introduced. Because the 
calculated density is renewed at each HF iteration, it is impractical to store the matrix 
elements of the density-dependent interaction. We here assume the usual zero-range form 
for the density-dependent interaction. The contribution of this component to the s.p. 
Hamiltonian is evaluated via the local densities. 

In reproducing the saturation, it could be an alternative way to introduce the momen- 
tum dependence in the central force, which satisfies the translational invariance, instead 
of the density dependence. Although we do not consider that possibility in this paper 
except for the case of the Skyrme interaction, the momentum-dependent two-body in- 
teraction will be handled in a similar manner to the momentum-independent interaction 
discussed below. 

In addition to the saturation properties which are relevant to the central force, the 
LS splitting is significant in the atomic nuclei. This suggests necessity of the LS and/or 
the tensor forces, though true origin of the LS splitting is still under discussion |T^. In 
most mean-field calculations so far, the zero-range LS force was assumed p|, |^. The 
zero-range tensor force was sometimes taken into account so as to cancel a certain term 



of the LS current ||T5[. Finite-range LS and tensor forces are also considered, as well as 
the zero-range ones. 

We thus consider the effective interaction in the following form. 



Vr. = T.it7+ttP'^-t'',Pr-tfP.Pr)f^{n2) 



V 



LS 



12 



E('^r^TE + t;;^^PT0)/f (ri2) Li2 ■ (Si + S2) 






^TN ^ 5:(t™-PT,+t™opTo)/r(n2)rL5i2 



DD 
^2 



t3(l + X3P.)[p(rir(5(ri2), (10) 



where /^ represents an appropriate function, /i stands for the parameter attached to the 
function {e.g. the range of the interaction), and t^ the coefficient. As examples of /^(r'12), 
the delta, the Gauss and the Yukawa forms will be considered. The relative coordinate is 
given by ri2 = ri — r2 and ri2 = |ri2|. Correspondingly, we define pi2 = (pi — P2)/2. Pa 
(Pr) denotes the spin (isospin) exchange operator. Pte (-Pro) is the projection operator 
on the triplet-even (triplet-odd) two-particle state, 

1 + P 1 — P 1 -I- P 1 + P 

2222 ^ ^ 

Similarly, the projection operators on the singlet states are defined by 

1 — P 1 + P 1 — P 1— P 

PsE = ^-^^-^, Pso = ^-^^-^. (12) 

L12 is the relative orbital angular momentum, 

Li2 = ri2xpi2, (13) 



Si, S2 are the nucleoli spin operators, and 5*12 is the tensor operator, 

Si2 = 4 [3(si • ri2) (S2 ■ ri2) - Si ■ Sa] . (14) 

The nucleon density is denoted by p(r). 

The Skyrme and the Gogny interactions are included in the present category of effec- 
tive interactions. In the case of the Skyrme interaction, we use fp{ri2) = S{ri2) for the 
momentum-independent central force. Discussions on the momentum-dependent and LS 
terms will be given in Appendices ^ and ^. In the case of the Gogny interaction, we set 
fuii^u) = exp[— (/iri2)^]. The LS force has the same form as in the Skyrme interaction. 
In both interactions, f ™ = is assumed, except for the counter-term to a part of the 
LS current. In some of the recent parameterization of the Skyrme interaction |T6[ the LS 



contribution is expressed only in the density-functional form, without explicit correspon- 
dence to the two-body interactions. They are not expressed in the form of Eq. ([TD|) and 
will not be considered in this paper. 

4 Calculation of two-body interaction matrix elements 



c 

125 



We now discuss how to compute the matrix elements of the two-body interactions, v 
v\2 and f ™. Their contribution to h^^^ in Eq. (Q) is given by 

E (Nn'.f) ,, . ^it l^U ^ ^S ^9 (i^^J^ arf)J\iv?2+V^!+vl2'')mJ, P'nV)A 

(15) 
where (Nnij) denotes the occupation number of the s.p. orbit. The expression (aij) 
represents the basis \'^aij) of Eq. (|l|) in Section ^j. Though not shown explicitly, the 
proton- neutron degrees-of-freedom should be considered in Eq. (|15[), in practice. As has 
been discussed in the preceding section, the contribution of vf^ to the s.p. Hamiltonian is 
calculated through the proton and neutron density distributions at each iterative process. 



as in Ref. [12 



Let us take the central force f^2 ^s an example. The anti-symmetrized two-body 
matrix element in Eq. (|T5|) is obtained from the non-anti-symmetrized ones, 

{iflJ2)J\v?2\ijlJ2)J)A = {ij[j2)J\v?2\UlJ2)J) " {{j[j2)J\v?2\ij2jl)J) , (16) 

where | )a (| )) denotes anti-symmetrized (non-anti-symmetrized) state vector. Without 
confusion, the symbol j is regarded as an abbreviation of (aij) and the proton- neutron 
degrees-of-freedom. As shown in Eq. (|T3|), we only need the anti-symmetrized matrix 
elements with ii = i[, £2 = i'2, ji = ][ and J2 = J2 fo'^ the spherical HP calculations. 
However, we here discuss how to evaluate the matrix elements in more general manner, 
as will be useful in the case that the symmetry is broken. Inserting Eq. (0), we have 

UJ'2)J\V'{2\{JIJ2)J) = Y.Uj2)J\(t^ + tlP. - f^Pr - tf P.Pr) f^ {ru)\{jlj2)J) ■ (17) 



Without loss of generality the spatial function f^{ri2) is assumed to be common for 
the Wigner, Bartlett, Heisenberg and Majorana terms. The difference among them are 
in the spin and isospin operators, for which we use the notation O^ (= 1 or Po-) and 
Or (= 1 or P-r). By converting the jj-coupling to the LS-coupling, each term of the 
interaction in Eq. ([T7|) is written as 



{{3[32)J\f'^{ri2)O^OM3,3,)J) 



Y,{2L + 1)(2S+ 1) J(2ji + l)(2j2 + 1)(2j1 + l)(2j^ + 1^ 
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J 



x{{l[QL\f;:{ru)\{hh)L){0^)s{Or), (18) 

where {Oa)s ((C^r)) denotes the expectation value of O^ {Or) of the spin (isospin) part. 
The spatial part in the right-hand side is defined by 

((^;^2)^l/,'^(n2)|(/i/2)i^) = /rfVirf\2i?,;(ri)i?,^(r2){[F(^'^)(r0r(^^)(r2)]i?}* 

x/,^(ri2)%(ri)P,,(r2)[r(^^)(fi)F(^^)(f2)]£^ (19) 

The spatial matrix element (^) can straightforwardly be calculated for simple forms 
of the interaction such as the delta form. However, we here intend to handle various 
types of interactions, including the Yukawa form. For this purpose, we utilize the Fourier 
transform of /^(ri2), as was exploited in Ref. [|^ for the HO bases. 

The Fourier transformation of f^{r 12) gives 

^(A;) = |rfVi2^(ri2)e-^^-^-. (20) 

By inverting this transformation, we obtain 



Mru) = J^sJ d'kf,{k)e^'^-^^ = ^^ / rf=^A;/,(A;)e^'^-^e-'^-^ . (21) 



Substituting Eq. (^) into Eq. ([19|), we find the ri and r2 integrals are separated at the 
expense of the k integration. The angular integration is implemented by using 



gik-r 



47r J2 ^^(2A + l)3x{kr) Y^^^ (k) • F^ (r) , (22) 

A 

where j\{x) denotes the spherical Bessel function, deriving 



((/l/^)L|/^C(ri2)|(/i/2)i^) = Ev(2^i + l)(2^2 + l)(^i0A0|£;0)(£20A0|f20) 

A 

xj^kHkf^{k)xf\k)l^^\k). (23) 
Here XJ and X2 are defined as 

■lf\k) = j^ r'drjx{kr)R^^,^,Xr)R.'A&) . (24) 

7 



with the subscript i (= 1,2) corresponds to the nucleon index, which actually represents 
{cidiji, Ci'i^iJi)- Since the radial part of the present basis Raij{r), given in Eq. (|l|), has the 
Gaussian form, X^'^^k) is calculated analytically. 



2:(°)(^) = 0,a'^' 



k \ ,(A+i) 



rhry/y' 



^ V- r 2 



k 



rirvry' 



exp 



k 



rXiCVfy' 



(25) 



where L^^\x) is the associated Laguerre polynomial and 



2 2 \-Pa+Pa' Y 



+ e -\ 



+ Pq + Pa' + 1 



Ca^,a'. 



^(2£ + l)!!(2f + 1)!! 



Z/„/ \ ^ I-Pq-Pq' / 2z/aZ/Q/ 



Vr, 



2 rPQ-|-Pc«/ 



.^^ + ^a'. 



r\jn/rv' 



V4+^- 



(26) 

(27) 



We thus obtain the following expression for the spatial part of the matrix element. 



m'2)L\fAri2)\{hh)L) 



5:v/(2£i + l)(2£2 + l)(40A0|fi0)(£20A0|f20)Ca,^,,„;,;C 

A 



Q!2«2,a24 



X / edkK(k) 



U7^Vl;^:|l 



_ i^aia' 



xL 



(A+i) 

12+1' -\ 

^ +Pc2+Pc' 



i^n^nL 



^ "^020^2 



2^ +Pci+P^/^ 

k \ 



'•^rvi rv' 



aia' 



exp 



^ "^QlCt'j^ 



^rbQjo^ 



(2^ 



Whereas it is not easy in general to evaluate numerically the multi-dimensional inte- 
grals to a high precision, we have only one-dimensional k integral in Eq. (pSl). Moreover, 
even this k integral is analytically carried out for the typical interaction forms. Recall 
that the associated Laguerre polynomial is defined as 



Tia + n + l) 



x^ 



" ^ ^ ^,T{a + q + l){n-q)\ q\ ' 



(29) 



The k integral in Eq. (^) turns out to be the sum of the integrals with the form 



dA;A;2"+2e-(fc/^)7^(^)^ 



(30) 



where n is a certain integer and k^ = {^/i^a^^a' + V'^a2a' ) ^^ ■^*^^ ^^^ zero-range interaction 
such as the momentum-independent term of the Skyrme interaction, fslr^) = S{ri2) leads 



to fsik) = 1. Here we substitute S for the suffix ^ to show the function form exphcitly. 
The integral of Eq. ( ^0]) therefore reduces to 

' dk fc^^+^g-e^/^)^ = (2^+1)" A -2n+3 _ (31) 

2"+2 ^ ^ 

For the Gaussian form such as the Gogny interaction, f^{ri2) = e~^^'^'^^^ derives f^{k) = 
{^ / nY e~'^^/'^^'^ . Thus the integral of Eq. ( PP] ) for the Gauss interaction is 



\ fl J Jo \ fi J 2"+2 \^k2 4^2 

In both cases the k integral of (pO|) yields an analytic function. For the Yukawa interaction, 
f^{ru) = e-^'^^^ / iiri2 leads to f^k) = 4:'k/i^{^i'^ + P). The integration of (|0D is still 
written in a compact form, by using the error function, 

/i JO /i + K 

= ^-^(-u'r ly(2r - m (-^\ - ^ e^^)^ Erfcf^ 
where 






Erfc(x) = / e-' dz . (34) 



As is shown in Appendix 0, the momentum- dependent interaction, such as contained 
in the Skyrme interaction, can be handled in an analogous manner. The treatment of the 
LS and the tensor forces are discussed in Appendices |B| and p. By the present technique we 
can deal with various interactions, either zero-range or finite-range, in a unified manner. 

In coding a computer program, we should prepare a subprogram for the integration 
of Eq. (pOD . This integral of Eq. ( |30D is the only part dependent on the interaction form. 
Therefore various interaction form can be handled just by substituting the subprogram. 
Moreover, it is unnecessary to carry out numerical integration in calculating the inter- 
action matrix elements, for the delta, the Gauss and the Yukawa interactions. Even 
for a more complicated form of the interaction, numerical integration is only needed for 
Eq. (pOD, as far as its Fourier transform /^(fc) is known. 

The above technique is also applicable to the Coulomb interaction. The interaction 
form of f{ri2) = l/ri2 yields f{k) = ATc/k"^, and the k integral of Eq. ( pO| ) becomes 

An / dk e^e-^^/~^^ = {In - l)!!7r2 . (35) 







This is immediately obtained from the yU — > limit in the Yukawa interaction. Although 
the Coulomb exchange energy among the protons was often approximated |ll5| as 



^er--^e^@'/Mr)]SrfV, (36) 



this approximation can be lifted in the present approach. It should be mentioned that an 
alternative method for exact treatment of the Coulomb energy was proposed recently |T^ , 
where the Coulomb force is transformed into an integration of Gaussians. 

Once storing the interaction matrix elements and having an estimate of c^^^} i^ Eq. @, 
we obtain the s.p. Hamiltonian from Eqs. (|^,|1^), and via the densities for the contribution 
of vf^. Solving the HF equation by iteration, we can implement the HF calculation. 

5 Numerical tests 

In this section we shall demonstrate the present method via the HF calculations for the 
oxygen isotopes, using both zero- and finite-range interactions. The HF calculations are 
carried out with maintaining the spherical symmetry and the parity conservation. If the 
valence orbit is partially occupied, the contribution of the orbit to the s.p. Hamiltonian 
is averaged over the magnetic quantum numbers m, as shown in Eq. ([151) . We first 
investigate the characters of the bases, by taking the Skyrme interaction with the SLy4 



parameter-set |]T9[. Secondly the application to a finite-range interaction, for which the 
Gogny DIS ^ is used, will be shown. In both cases the Coulomb interaction is exactly 
treated, as mentioned above. The center-of-mass energies are corrected approximately, 
by taking only the one-body kinetic part into account. 

It takes longer time to implement HF calculations in heavy nuclei than in light nuclei. 
Though we restrict our application to the oxygen isotopes in this paper, it is sufficiently 
practical to apply the present method to the Pb isotopes, as will be shown in Ref. ||^. 

5.1 Selection of single-particle basis-sets 

We use several sorts of single-particle basis-sets, composed of the bases having the form 
of Eq. (|I|). Each of the bases is characterized by the index a, which actually corresponds 
to the parameters Pa and z/q,. In practical calculations, we restrict the values of Pa and 
z/q, to a certain extent; otherwise there are too many possibilities. As mentioned earlier, 
we can take both the HO-equivalent basis- set and the KG basis- set, by choosing the 
parameters appropriately. As well as these sets, a hybrid basis-set will be tested, in 
numerical calculations shown in the subsequent subsections. 

The basis-set equivalent to the HO one is obtained from Eq. (|l|), by posing 

Pa = a-1, Ua = Uoj, {a = l,2,---,K) (37) 

where 

^ (38) 



^ \ Mu 

In the numerical calculations shown below, we do not consider the nucleus dependence of 
the Vi^ parameter, taking Tiuj = 41.2 x 24~^/^ MeV, for the sake of simplicity. After the 
Gram-Schmidt orthogonalization, the basis index a corresponds to the number of nodes 
n in the HO bases. In the usual calculations by the HO bases, the truncation is made in 
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terms of A'^sh = 2n + L However, we here fix the number of the s.p. bases K irrespective 
of i and j, to be fair with the case of the KG set mentioned below. This indicates the 
truncation according to ra, rather than by A'^sh- In the above set of (0), the maximum 
value of n corresponds to i^ — 1. 

As stated in Section ^ the bases of Eq. (|1]) are not orthogonal and care must be taken 
so that the norm should not be too small after the orthogonalization. If we adopt the HO- 
equivalent set of Eq. (p7|), the norm of the i^-th basis appreciably decreases for growing 
K. Numerical instability seems to occur for i^ > 10 in computations with the double 
precision. Hence we always restrict ourselves to K = 7 when we use the HO-equivalent 
basis-set. 

The KG basis-set is obtained by 

Pa = 0, z/„ = i/i&"-i. (a = l,2,---,ir) (39) 

We use the same vi, b and K for all i and j. If the common ratio b is close to unity, 
the overlap between the a-th and the {a + l)-th bases is large. Then the norm after the 
orthogonalization becomes vanishingly small, which may lead to numerical instability. On 
the other hand, if we adopt the larger value of b, it is the more difficult to reproduce the 
wave-functions accurately. For instance, in order to represent the smooth exponential 
decrease of the density by a superposition of the Gaussians, b should not be very large. 
In practice, the density distribution shows bumpy structure for b > 1.35. In the following 
calculations we fix 6 = 1.33, so as for the exponential asymptotics to be reproduced in 
an effective manner. For the range parameter Uq., we take one of them to be equal to i/^ 
in Eq. (^8]). The z/i value is determined accordingly; for example, z/i = z/^ 6~^ if we set 

In the HF calculations we always confirm the convergence for iteration. However, it is 
not easy in most cases to pursue the convergence for increasing K. This is also true for 
the KG basis-set. The KG basis-set is characterized by three independent parameters; the 
shortest range z/i, the longest range ux and the common ratio b. Correspondingly, there 
are three courses to increase the number of bases K. One is to add the bases i'k+i and 
so forth, which have longer ranges than uk, with fixed z/i and b. The longer-range bases 
might be important to reproduce the wave-functions in the asymptotic region, particularly 
for the drip-line nuclei. Another is to shorten ui, keeping the longest range and b. The 
shortest range z/i is primarily relevant to the wave-functions deeply inside the nucleus. 
The other is to take smaller b. This could be significant to accurate description of the 
wave-functions in any region. In order to attain the full convergence in the HF calculation, 
all of the three courses should be tested. 

While the KG set has an advantage in describing the wave-functions in the asymptotic 
region, it depends on the parameters how well the wave-functions in the surface region is 
reproduced. For the better description of the surface region by the KG set, we usually 
need the smaller b. An alternative way may be given by a hybridization of the KG set 
and a small number of the HO-type bases. In addition to the HO-equivalent and KG 
basis-sets, we also test the hybrid basis-set. Among the bases for each i and j, {K — 1) 
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bases are taken to be the KG ones, and for the last basis we use pa = ^', 

= 0, u^ = u,b»-\ (a = l,2,---,i^-l) 




' " ' ' ^ ' ' ' ^ (40) 

The parameters b and z/^^ are taken to be the same as in the KG and HO bases mentioned 
above. 

5.2 Case of zero-range interaction 

We apply the present method to the even-iV oxygen isotopes. Using the Skyrme SLy4 
interaction, we first compare results from the HO (the HO-equivalent, in practice), the 
KG and the hybrid basis-sets for a fixed value oi K; K = 7. For the KG and hybrid 
sets, we take ui = u^h'^, leading to the longest range vk — 5.7 fm for the KG set and 
uk-i — 4.3 fm for the hybrid set. 

The variational character of the HF theory is available in comparing results among 
different basis-sets. As the total energy is lower, it is closer to the true HF energy, in 
principle. The total HF energies calculated with the HO, the KG and the hybrid basis-sets 
are shown in Fig. |l|, where the energy differences from the lowest one is plotted for each 
nucleus. In comparison with the HO set, the KG set gives higher energies for the A < 22 
oxygen isotopes, while in A > 24, where the neutron lsi/2 orbit is occupied, the KG set 
gives lower energies than the HO set. This is ascribed to the broad radial distribution 
of the lsi/2 orbit, which is hardly reproduced by the HO basis-set. The hybrid basis-set 
works very well in the whole region. The energies are close between the HO and the 
hybrid sets in ^^~^^0, having the differences less than 0.01 MeV, and the hybrid set yields 
sizably lower energies than the KG set for all of the calculated oxygen isotopes. Thus the 
hybrid basis-set of Eq. (^OD is adaptable both to stable and unstable nuclei. 

In Fig. ^, the density distribution is compared among the three basis-sets, for ^^O, 
^^O and ^^O. The density distribution in r > 6fm tends to obey to the exponential 
asymptotics. Obviously, the density by the HO basis-set does not distribute sufficiently 
as A increases, unable to reproduce the asymptotics for ^^O and ^^O. The densities in 
r > 6fm calculated with the HO set behave quite analogously among the three nuclei, 
suggesting that they originate in the character of the bases and are not physical. Therefore 
the HO set is practically incapable of reproducing the asymptotics. On the contrary, the 
KG and the hybrid basis-sets reproduce the exponential asymptotics rather well. Although 
the KG set does not give the exact exponential asymptotics, it is possible to approximate 
the asymptotics by the KG set in an effective sense. The same holds for the hybrid set. If 
the density becomes extremely low, it is difficult to be reproduced by the KG set. For this 
reason, the KG set shows fictitious behavior of the density for r > 9fm in ^^O, although 
the hybrid set yields rather smooth decrease. For ^^O and ^^O, the densities obtained by 
the KG set distribute more broadly than those by the hybrid set. This seems to caused by 
the difference in the longest range of the bases, on which it depends how slowly decreasing 
density can be described. Because we set the longest range to be u^^ It" in the KG set while 
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Ui^ IP' in the hybrid one, the KG set can reproduce broader distribution that the hybrid 
set. 

Figure ^ shows the high adaptabihty of the KG basis-set, particularly for the wave- 
functions in the asymptotic region. The energy-dependent asymptotics are reproduced 
reasonably well by a small number of bases. Unlike the HO basis, an individual basis in 
the KG set will not be a good first approximation of the nuclear s.p. wave-function. Nev- 
ertheless, when a certain number of bases are superposed, the KG set acquires remarkable 
flexibility in describing wave-functions. In the nuclear wave-functions, the exponent of the 
asymptotic form is energy dependent, as is viewed in the nucleus dependence in Fig. |^. 
The KG set (and therefore the hybrid set) well approximates the exponentially damping 
wave-functions with a moderate number of bases, whatever the separation energy is. 

We next increase K, the number of the bases for each £ and j. The HF calculation 
is carried out using the KG and the hybrid basis-sets for i^ = 10 and 15. In the case 
oi K = 10, we take vi = u^ b~^ both for the KG and the hybrid sets, giving the longest 
range uk — 10 fm for the KG and i'k-i — 7.5 fm for the hybrid set. In the calculation 
with K = 15, we take ui = v^^h~^ for the KG set and z/^ 6~^ for the hybrid set, both 
having the longest range of 31 fm. In Fig. ^, the HF energies and the rms matter radii are 
plotted as a function of K, for ^^O, ^^O and ^^O. The density distributions obtained by the 
hybrid set of i^ = 15 are already shown in Fig. ^. For the matter radii, the center-of-mass 
correction is neglected, corresponding to the density distributions depicted in Fig. ||, in 
order to view the properties of the basis-sets directly. Owing to the variational character, 
the HF energies become lower as K increases. If we adopt the KG set, the HF energies 
decrease slowly for increasing K. On the contrary, the HF energies by the hybrid basis-set 
are stable between i^ = 10 and 15, where the biggest decrease among ^^-^sq jg ]2ierely 
0.003 MeV. This sort of stability is also viewed in the rms matter radii. Although the 
radii tend to be underestimated by the hybrid set with K = 7, their difference between 
the K = 10 and 15 cases are negligibly small; less than 1.5 x 10~^fm for all the 1'^~28q 
nuclei. By comparing with the K = 15 hybrid basis-set, we view that the HO set provides 
reasonable values of the HF energies and the matter radii in the stable nucleus ^^O, while 
it is not satisfactory in the unstable ones such as ^^O and ^'^O. 

Though one may think that the HF energies are convergent in the K = 15 result using 
the hybrid basis-set, it does not imply the full convergence. The HF energies further 
decrease, if we use smaller b. Indeed, the HF energy becomes lower by about 0.01 MeV 
for a few of the oxygen isotopes, if we use b = 1.25. As in the calculations so far, it is quite 
a laborious work to attain the full convergence in the HF calculation, and is beyond the 
scope of this article. In the same regard, the HF energies shown in Figs. |l| and |^ do not 
immediately mean a drawback of the KG set itself. It is only indicated that the KG set 
with 6 = 1.33 is not sufficient to describe the nuclear wave-functions in the surface or the 
interior region. If we use smaller b than the present value, the wave-functions around the 
surface may be reproduced more accurately, though it requires larger number of bases. It 
will be fair to say that the convergence is accelerated by using the hybrid set, compared 
with the case of the KG set. 

In the present approach, the most time-consuming part in the numerical calculation 
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is the computation of the two-body interaction matrix elements. It takes about 430 sec of 
CPU time on HITAC SR8000 to calculate all the necessary matrix elements up to the sd- 
shell when K = 7, although the program has not yet been tuned. Moreover, the CPU time 
for computing the two-body matrix elements is almost proportional to K^. On the other 
hand, the HF iteration need about 30 sec for each nucleus, and has dependence weaker 
than K^. Under this situation, a great advantage of the KG basis-set is that it does not 
include parameters specific to mass number or nuclide. Hence the same basis-set can be 
used for a number of nuclei. With this advantage of the KG set, we have to calculate the 
two-body interaction matrix elements only once and do not have to recalculate them in 
systematic studies, and thereby we can save the computation time. 

5.3 Case of finite-range interaction 

The present method of the HF calculation can be used for finite-range interactions. We 
demonstrate it via the calculation with the Gogny DIS interaction. 

There is a problem in the Gogny DIS force when it is applied to the mean-field calcu- 
lations using the KG basis-set |^. For the pure neutron matter with the DIS force, the 
energy per nucleon diverges with the negative sign at the high density limit. Originating 
in this defect, the HF energy goes to negative infinity in the finite nuclei, when all the 
neutrons gather in the vicinity of the origin {i.e. the center-of-mass) without overlap 
of the proton distribution. For the /3-stable nuclei, this unphysical configuration is well 
separated from the normal HF solution, i. e. an energy minimum satisfying the saturation 
properties, and the normal solution is stable enough to be obtained in the numerical cal- 
culations. However, it is not the case for the highly neutron-rich nuclei. Even if the initial 
configuration is in the physical domain, the tunneling to the unphysical configuration 
takes place before convergence. In order to circumvent the tunneling, we need a certain 
cut-off of the high momentum components. In the previous studies [^, a sort of cut-off 
was implicitly introduced by adopting a limited number of the HO bases. On the other 
hand, we have to be cautious when we use the KG set. The wave-function of ^^O collapses 
via the tunneling when bases having Ua < lim. are included. Alternative to cutting off 
the high momentum components, a way to avoid this problem is to modify the interaction 
parameters, say, X3 in vf-P. This possibility will be explored in a forthcoming paper [0. 

In the numerical calculations, we use the HO, the KG and the hybrid basis-sets of 
(P7|), (P^) and (^D|). For the HO set, we take K = 7 ioi each £ and j, assuming the same 
value of z/(^ as in the preceding subsection. This corresponds to the A^sh < 13 truncation 
except for the (i-orbits, for which the Agh = 14 bases are included, and this basis-set is 
similar to that used in most mean-field calculations with the Gogny interaction so far. 
For the KG and hybrid sets, we use ui = v^^h'"^ and h = 1.33, i.e. vi = 1.36 fm, and 
K = 12. 

The density distributions of ^^O, ^'^O and ^^O are depicted in Fig. ^. The KG as 
well as the hybrid bases yield the reasonable asymptotics for unstable nuclei, whereas 
the HO set does not. The variational character is not fully available until establishing a 
rigorous cut-off scheme. However, it is still informative to compare the HF energies, as 
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shown in Table |1|. It is very likely that the energy of ^^O by the HO set is hampered 
by the ill asymptotic behavior, being higher than the energy of the KG set. As in the 
SLy4 case, the hybrid set gives the lowest energies in all of ^^O, ^^O and ^^O. Thus the 
wave-functions in the asymptotic region, which have a sizable contribution to the physical 
quantities in highly neutron-rich nuclei, might not be described properly in the previous 
Gogny mean-field calculations by the HO bases, even though the i/^ parameter is better 
tuned than in the present calculation. 

We next compare the HF results of the Gogny DIS interaction with those of the 
Skyrme SLy4 interaction. In Fig. ^, the neutron s.p. energies of the S(i-shell orbits are 
shown for the oxygen isotopes. For ^'^"^''O, the 0^3/2 orbit is unbound in the DIS results, 
and hence its energies are not presented. The s.p. energies obtained from the SLy4 and 
the DIS interactions are relatively close to each other. A notable difference is found in the 
behavior of e„(lsi/2) around ^^O; the DIS interaction gives a kink at ^"^O, while the SLy4 
does not. When the N (neutron number) dependence of the s.p. energies in 16-22Q jg ^^ 
effect of the occupation of 0^5/2, changes in the s.p. energies from ^^O to ^"^O and from 
^^O to ^^O are connected to the lsi/2 and 0^3/2 occupation. It was argued based on the 
recent experimental data p^ that A^ = 16 becomes a magic number in the neutron-rich 
region. The behavior of the s.p. energies around ^^O could be relevant to the magicity of 
A^ = 16. The kink in the DIS result gives rise to the relatively large gap between eri,(lsi/2) 
and 6^(0^3/2) at ^"^O. This would make the ^^O core stiffer than in the SLy4 interaction. 
We have confirmed that the kink viewed in the Gogny DIS result does not emerge in the 
results of other popular parameter-sets of the Skyrme interaction, as well as in the SLy4 
result. Still it is not clear whether or not the range of the interaction plays a role in this 
kink. It is also commented that the kink at ^^O in the DIS interaction is not apparent 
when we use the HO basis-set, probably because of the wrong asymptotics. 

The two-neutron separation energy S2n is calculated as difference of the binding en- 
ergies between the neighboring isotopes. The calculated values of 5*2™ by the hybrid 
basis-set [K = 15 for SLy4 and A' = 12 for DIS) are shown and compared with the 
measured ones [^ in Fig. |^. The 5*2^ values are not very different between the SLy4 and 



DIS interactions for ^^ ^^O. It has been confirmed experimentally that ^^O and ^^O are 
unbound 123, E3|. This indicates 5'2„ < for ^^O. This feature is not reproduced in the 



SLy4 result. In the DIS interactions, we have slightly positive S'2„. However, this seems 
to depend somewhat on the details of the numerical set-ups; for example, the treatment 
of the center-of-mass correction. We just state that, in respect to 5'2„, at ^^O, the DIS 
interaction gives preferable result to the SLy4 interaction. 

In the present calculations we do not take into account sufficiently the collectivity 
due to the pairing correlations. If we rely on the shell closure at A^ = 16, the pairing 
correlations will hardly change the energy of ^^O, while they lower the energies of ^^O and 
^^O to a certain extent. Thus the pairing effects are expected to give lower 5*2™ at ^'^O and 
higher 6*2,1 at ^^O than the present HF values, if we perform a Hartree-Fock-Bogolyubov 
calculation. 

Concerning the computation time, it is noted that the interaction dependence is weak 
in the present method. The CPU time for the Gogny interaction is almost the same as in 
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the Skyrine interaction. 

6 Summary and outlook 

We have developed a new method of the Hartree-Fock calculations. This method has 
advantages in reproducing the slowly decreasing density distributions in unstable nuclei 
by a proper selection of bases, and in treating various interactions including finite-range 
ones. 

The key point of the method is adoption of the Gaussian bases shown in Eq. (|l]). 
This covers the Kamimura- Gauss (KG) basis-set, as well as the basis-set equivalent to the 
harmonic-oscillator (HO) one, and may open wider variety. We have also discussed a way 
to calculate two-body interaction matrix elements, by applying the Fourier transforma- 
tion. This is particularly suitable to the Gaussian bases because the numerical integration 
can be avoided to a great extent. Owing to this treatment of the interaction, we can easily 
switch from an effective interaction to another. 

The present method has numerically been tested by using the Skyrme SLy4 and the 
Gogny DIS forces, as representatives of the zero- and the finite-range interactions. The 
calculations with the HO basis-set and with the KG basis-set are compared. It has been 
confirmed that the KG set efficiently describes wave-functions in the asymptotic region for 
the neutron-rich nucleus such as ^'^O and ^^O. When we adopt the KG set, we do not have 
to change the bases from nucleus to nucleus, since they do not contain nucleus-dependent 
parameters like hu. Hence the KG basis-set is expected to be powerful for systematic 
calculations. We have also shown a way to improve the convergence over the KG set. If 
we use a hybrid basis-set, in which a HO-type basis is added to the bases in the KG set, 
the HF energies often decrease substantially. The results on the s.p. energies and on S2n 
are also compared between the SLy4 and the DIS interactions, for the oxygen isotopes. 

The present method provides us with a useful tool to investigate structure of the 
unstable nuclei, particularly with finite-range interactions. It will also be interesting to 
reconsider the effective interaction within the mean-field approaches. We can deal with 
various finite-range interaction, not only for the central part, and even with the Yukawa 
form. A research project in this line is under way. 

While we have assumed the spherical symmetry in the discussions in this paper, it is 
straightforward to extend it to the deformed nuclei. The future plan includes the extension 
to the Hartree-Fock-Bogolyubov approach, and the combination with the complex-scaling 
method so as to handle the resonant single-particle orbits. 

The authors are grateful to K. Kato and H. Kurasawa for helpful discussions. This 
work is supported in part as Grant-in- Aid for Scientific Research (C), No. 13640263, by 
the Ministry of Education, Culture, Sports, Science and Technology, Japan. Numerical 
calculations are performed on HITAC SR8000 at Information Processing Center, Chiba 
University. 
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Appendices 

A Matrix elements of momentum-dependent part of 
Skyrme interaction 

The central part of the Skyrme interaction is parameterized as 

vg = to(l + xoP^)6{ru) + -ti(l + XiP^)[6{ru)Pu + P?2'^(ri2)] 

+ t2(l + a;2Pa)pi2-5(ri2)pi2. (41) 

As has been discussed in Section |[ the two-body matrix elements of the to term can be 
treated in a unified way with the finite-range interactions. We here show that, though 
the ti and ^2 terms depend on the relative momentum pi2, they can also be treated in a 
similar manner by using the Fourier transformation. 
We first recall the identity 



[pi2,[pi2,5(ri2)]] = -(Vf2^(ri2)j , (42) 

where V12 = (Vi — V2)/2. By using the Fourier transform of the delta function, we obtain 

[Pi2, [Pi2, 5(ri2)]] = [P?25(ri2) + 5(ri2)p?2] - 2pi2 • 5(ri2)pi2 = 7^ / A;2e^k.r,.^3^ _ (43) 

Owing to the delta function, [p^2'^(i^i2) + ^(ri2)Pi2] vanishes when it operates on the 
spatially odd two-particle states. Similarly, 2pi2 • S{ti2)Pi2 vanishes when acting on the 
spatially even states. Hence we can write 

[pi,5{ru) + 6{ru)pl2] = ^ / ^V^-^-d^A: ■ (Pse + Pte) , 

-2pi2 • 5(ri2)pi2 = j^lkV^^-'^^'d^k-iPso + PTo). (44) 

The projection operators Pse, Pte, Pso and Pto can be incorporated in the spin-isospin 
part, by using Eqs. (|TI]JT^. Then the spatial matrix elements of [p^2'5(''i2) + '^(ri2)Pi2] 
and — 2pi2 ■ 5(ri2)pi2 are both evaluated by setting f^i{k) = k"^ in Eq. (pH]). 

B Matrix elements of LS interaction 

The LS interaction f|"| in Eq. ([T0|) can be handled in a similar manner to the central 
force. We here consider the non-anti-symmetrized matrix elements of the LS force, 

{{j[f2)M^\ijlj2)J) = E(Ui2)^l(^r^TE + ti^°PT0)/f (ri2) U2 " (Sl + S^W^J^U) • 

(45) 
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The LS force operates only on the spin-triplet two-particle states. As in Appendix 0, 
we separate Pte and Pto; denoting the projection operators by 0^r{= Pte or Pro) and 
their expectation values by {Oo-t)s=i- It is noted that /^^(ri2) should not be the delta 
function, since (5(ri2)Li2 = 0. 

From the definition of Eq. (|l3), L12 is rewritten as 



L12 = -(^1 + 4 - ri X p2 - r2 X pi) , 



(46) 



where £1 = ri x pi and £2 = 1*2 x p2. Since the (£1 + £2) operator does not change the 
spatial part of the wave- functions, the matrix elements regarding {£1 + £2) is handled in 
an analogous way to the central force. 



{{flJ2)J\f,iri2) Ah +£2) ■ (Sl + S2)a.|(jU2)J) 



E 3(2/^ + 1) 

L 



J(J + l)-L(L + l)-2 



(2ji + l)(2i2 + l)(2j( + l)(2j^ 



X 



h 2 31 
h \ 32 
L 1 J 



2 Jl 



^1 ^ - ■ 

1'2 \ 32 \m2)L\lTiri2)m2)L){0.r) 

L 1 J 



5=1 



(47) 



The {{l[l2)L\fj^^{ri2)\{hl2)L) matrix elements are given in Eq. (p8[) , except that f^{k) is 
replaced by /^^(A;), the Fourier transform of fh^{ri2). 

For the part including (ri x P2), we separate the spatial and the spin parts again, 
having 

{{3[32)J\fiHn2) ^(ri X P2) • (Si + S2)a.|(jij2) J) 



E 3V6(2L + 1)(2L' + l)(2ji + l)(2j2 + l)(2jl + l)(2j^ + 1) WiLJll;lL') 

L,L' 



h \ 3\ 
I2 \ 32 

L 1 J 



V - i' 
n 2 Ji 



1 



1'2 \ f2 }{il'il2)L'\\fi'in2)-{r,xp. 
V 1 J 



\m\ 



ll2)L) 



x{C)s=i, (48) 



As has been shown for the central part in Section ^, /^^(ri2) contains the angular part 
[y'^^^(fi) ■ Y^^\r2)]- Combining it with the (ri x P2) operator, we obtain 



FW(fi).FW(f2)] -(rixp2)« 



\A+1, 



'2A + 1 



rW(ri)rW( 



r2 



(0) 



)'-"' E /'^^'^r^^'^ A010|A,0)iy(llA,A;lA2) 



Al,A2 



xr.JF(^i)(ri) [FW(f2)V 



7(1) 



' 1 V 2 



(A2)l« 



(1) 



• (49) 
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The Fourier transformation of Eq. pO| ) and the integration of the angular part yields 



{(''i4)i'll/f('-i2)5(r,xp2)<')||((,fc)L) 



2 (_)>.,2A, + W5PAi±i)P^±i)P£±i) (A 1 0|A. 0) ,y (1 1 A. A; 1 A,) 

Ai,A2 * ^ 




Wi 



k^dk fj:^{k)^2h + 1 (/i Ai 0|/; 0)ll'{k) 

X <{ ^(^2 + l)(2/2 + 3) (/2 + IO A0|/^0) W^(/2 1 /2 A; /2 + I A2)X^'+^(A;) 

-^12(212 - 1) ik-lO X0\l20)W{l2ll'2X;l2-l X2)4'''\k)\ , (50) 



where 



xf^(A:) = 


= / r^drrjx{kr)Rj>^r) 

^ 


%W, 




r\k) - 


= / r'^drjx{kr)Rj'{r) 
Jo 


\( d i\ 1 


■) 


r\k) - 


= / r'drjx{kr)Ry^{r) 
^ u 


\dr r J ' 


Xr) 



(51) 



The subscript i to X corresponds to the nucleon index. The Gaussian bases of Eq. (|l|) 
give 



dr 



d__t 
dr r ^ 

+ v 



Raej{r) 
Raijir) 



'2po 



r 



2r 



vt: 



R. 



aej\ 



'2i + 2p^ + l 2r\ ,, 
^jR^iji^)- 



(52) 



Hence the r integration in Eq. (51) is implemented in a similar manner to Eq. (pSQ . 
Because of the parity selection rule, the result is a product of a polynomial of k and a 
Gaussian, and finally the k integration has the form of Eq. (^Ol). Thus the k integration 
is carried out analytically for /^^(ri2) = e"^^*"^^^ , while it is expressed by using the error 
function for /^^(ri2) = e~'^^'^^ /fir 12. The matrix elements of the part coming from (r2 x pi) 



in Eq. (46) are obtained by interchanging the nucleon indices 1 and 2, with an appropriate 
phase factor. 

In the Skyrme and the Gogny interactions, the LS part is taken to be 



v^^ = 2iWo [pi2 X 5(ri2)pi2] ■ (si + S2) . 



(53) 



This type of the LS interaction is also treated in an analogous manner, by using the 
Fourier transformation. It is noted that this f}f force operates only on the S" = 1, T = 1 
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lun U^\ e-^'^^-)' 



(54) 



[i.e. spatially odd) channel, because of the (si + S2) operator and of S{ri2). Recall that 

r 

Substituting Eq. (^3|) into Eq. (|^), we obtain 

3 



LS 

h2 



2iWo lim 
-AWo 



5 



P12 X e-^'^^-^'pia 



[S1 + S2) 



lim il-e-('^'-i2)' 



/i— ►oo 



7r2 



L12 ■ (si + S2) 



(55) 



We shall take the fx ^ 00 limit after a certain algebra. By comparing Eq. (|55|) with 
Eq. (|TUp, we can identify 



mru) 



-4 lim -^e-(^''i2)' , 



whose Fourier transform is 



/LS(A;) = _4 1im/i2e-('^/2M) 



/J— >oo 



(56) 



(57) 



and t^»° = M/q- We here use the suffix 6" to stand for the spatial part of the LS force 
(^), for the reason clarified in the discussion below. By expanding the exponential factor 
and taking the /i -^ 00 limit except the non-vanishing terms, we have 



/LS(A;) = -4 hm fi')+Ii'. 



(5^ 



Although the first term in the right-hand side looks divergent, it is obvious that its 
inverse Fourier transform is proportional to 6{ri2), which turns out to vanish because 
6{ri2) L12 = 0. We now pose safely, for the LS force of Eq. (|53|), 



f^m = k''- 



(59) 



This is the same form as in the momentum-dependent term discussed in Appendix |^. 
Equation (|59D corresponds to the modification of Eq. (|56|) as 



fhHru) = -^ }im^f^' 



/^ 



-(Mri2)^ _ 



6{r 



12 



(60) 



Thus the LS force in the Skyrme and Gogny interactions is treated in a unified way with 
the finite-range LS forces. 
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C Matrix elements of tensor interaction 

We turn to the non-anti-syminetrized matrix elements of the tensor force, 

{{J[J2)J\VJ^\{J1J2)J) = E{ijlJ2)m7''PTE + tr'^PT0)/r(-12) r?2^12 | (jlj^) J) • (61) 

The tensor force operates only on the spin-triplet two-particle states. As in Appendix 0, 
we shall denote Pte or Pto by Oar- 

By separating the spatial degrees-of-freedom from the spin ones, the tensor operator 
of Eq. (jn^ is rewritten as 



Gtt r 



rl^Su = 8UI^ 



r^r(')(fi) + r2r(2)(f2)l -47rrir2 \Y^^\h)Y^^\r2 



(2) 



■Ji)Ji) 



The matrix element of the rfY^'^\ri) part is given by, after calculating the spin part. 



(2) 

(62) 



((j;j^)^l/r(-i2)-8J|:r?F(^)(fO- 



■bi 62 



(2) 



0.r\Ulj2)J) 



= Yl 12V5(2L + 1)(2L' + l)(2ji + l)(2j2 + 1)(2j1 + l)(2j^ + 1) H^(L J2 1; 1 L') 

L,L' 

X< /2 i J2 1'2^ \ f2 ((/'l/2)^1l/r(^12)\/Y^?^^'nfl)ll(^1^2)L) 

L/ L J \ \ L/ i •-'J 

x(ar)5=i. (63) 

By using the Fourier transform of /™(ri2) and integrating out the angular part, we obtain 
for the spatial matrix element, 

Qtt 



m2)L'\\C\n,)J-rtY^^>ih)\m2)L) 



E 

Ai 



'3(2Ai + l)(2L + l)(2L'+l) 



h h L 
(A020|AiO)<( Ai A 2 

n '2 ^ 



I' n\ 'r(2) 



X / k'dkf'''{k)J2h + l{h0Xi0\l[0)Il'>{k)-J2l2 + l{l20X0\l20)ir{k), (64) 



/' n\T(o)i 



where 



jf)(A;) = r r'^drr^jxikr)Rj>{r)RjXr) 
^ u 



(65) 



and X\ is defined in Eq. (^4]). The matrix element of the r2y*^^-'(f2) part in the expansion 
of Eq. (|6^ is obtained by an appropriate replacement between the nucleon indices. 
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After the algebra regarding the spin degrees-of-freedom, we have for the matrix ele- 
ments of the rir2[F*^^^(ri)F*^^^(f2)]*-^'' part, 



(2) 



(2) 



Y^ l2Jh{2L + 1)(2L' + l)(2ji + l)(2j2 + 1)(2j1 + l)(2j^ + 1) W{L J 21- 1 L') 



L,L' 



h \ Jl 



/' - i' 
n 2 Ji 



h \ 32){ 1'^ \ 3'2 ((/;/^)L'||/™(ri2)47rrir2 Y^'\h)Y'^'\h 



L 1 J 



L' 1 J 



(2) 



1/2)^) 



x(ar)5=i. (66) 



For the spatial matrix elements, we obtain 

((^'1^2)^'! 1/^(^2) 47rnr4F«(ri)r«(r2) 



(2) 



1/2)^) 



= {-Y J2 3^5(2Ai + l)(2A2 + l)(2L + l)(2L' + l)(A010|Ai0)(A010|A20) 

Ai,A2 

( h h L \ 

xVr(21AiA;lA2)< Ai A2 2 \ fc^rffc^M + 1 (/lO Ai 0|/;0)t{^^(A;) 

x^2/2 + l(/20A20|/;0)J?^(A;), (67) 
where Xi is given in Eq. (|5T|). 



In Ref. |T^, a zero-range tensor interaction was proposed to cancel out a certain term 
of the LS current. By separating it into the spatially even and odd channels, we have 



,TN 
h2 



■12 



AtlJifF- { [3(si ■ pi2)(s2 ■ P12) - (si • S2)Pi2] '^(ri 

+(5(ri2) 3(Si ■ Pl2)(S2 • P12) - (Sl • S2)Pi2] } ^TE 

-AtlJi^P {3(si ■ Pi2)5(ri2)(s2 ■ P12) + 3(s2 ■ Pi2)5(ri2)(si ■ P12) 

-2(Si • S2) [Pl2 ■ (5(ri2)pi2]} i^TO • 

In Ref. |T^ t™/° = |t™/^. By using the expression of Eq. (0), we obtain 



/^:^(ri2) 



lim 



4/1^ 



<^^rl2f 



TT2 



whose Fourier transform yields 



/<S"" (^) 



(6^ 



(69) 



(70) 



after eliminating the vanishing terms. Thus this type of the tensor force is also treated in 
a unified way with the finite-range tensor force. 
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Table 1: Hartree-Fock energies (MeV) in the DIS interaction for ^^O, ^^O and ^^O, using 
the HO, KG and hybrid basis-sets. 

nuchde HO KG hybrid 



1^0 -136.849 -136.720 -136.952 
2^0 -179.477 -179.445 -179.815 
2^0 -181.397 -181.501 -181.770 
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Figure 1: Comparison of the calculated HF energies in the SLy4 interaction among the 
HO, the KG and the hybrid basis-sets for the oxygen isotopes. The circles, pluses and 
crosses indicate the energies obtained from the HO, KG and hybrid sets, respectively. In 
all cases K = 7. In each nucleus, the energy differences are measured from the lowest HF 
energy of those by the three basis-sets. 
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Figure 2: Comparison of the density distribution in the SLy4 interaction among the HO, 
KG and hybrid basis-sets for ^^O, ^^O and ^^O. The dashed, dot-dashed and thick sohd 
hues are obtained from the HO, KG and hybrid sets with K = 7, respectively. As a 
reference, the densities obtained with the K = 15 hybrid set are shown by the thin sohd 
hne. 
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Figure 3: HF energies and rms matter radii of ^^O, ^^O and ^^O in the SLy4 interaction, 
as a function of K {i.e. number of the bases). The circles, pluses and crosses represent 
the results of the HO, KG and hybrid basis-sets, respectively. The dashed lines are drawn 
to guide eyes. 
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Figure 4: Comparison of the density distribution in the DIS interaction among the HO, 
KG and hybrid basis-sets for ^^O, ^^O and ^^O. The dashed, dot-dashed and sohd hues 
are obtained from the HO, KG and hybrid sets, respectively. 
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Figure 5: Neutron single-particle energies for the oxygen isotopes, obtained from the HF 
calculations with the SLy4 (open circles) and with the DIS (diamonds) interactions. The 
lines are drawn to guide eyes. 
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Figure 6: Two-neutron separation energies for the oxygen isotopes, obtained from the HF 
calculations with the SLy4 (open circles) and with the DIS (diamonds) interactions. The 



cross symbols represent the experimental data [^ 
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